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The formation of conducting channels of TuOr inside Ti02-based memristors is believed to be 
the origin for the change in electric resistivity of these devices. While the properties of the bulk 
materials are reasonably known, the interface between them has not been studied up to now mostly 
due to their different crystalline structures. In this work we present a way to match the interfaces 
between TiC >2 and Tu 07 and subsequently the band offset between these materials is obtained from 
density functional theory based calculations. The results show that while the valence band is located 
at the TL 4 O 7 , the conduction band is found at the TiC >2 structure, resulting into a type II interface. 
In this case, the Tu 07 would act as a donor to the Ti02 matrix. 


PACS numbers: 

I. INTRODUCTION 

The memristor is an electronic device predicted by L. 
Chua in 1971 [H H| and experimentally obtained a few 
years ago. [3] These devices are characterized by distinct 
resistance states which could be used to store informa¬ 
tion as they show very fast switching and high retention 
times. [4] 

There are reports of a wide range of materials (particu¬ 
larly metal oxides) which show this memristive effect. [E|- 
0 Although the switching mechanism, and the origin 
of the different resistance states has not been completely 
elucidated, it is generally believed that the presence of 
oxygen deficient phases plays a major role in the mem¬ 
ristive effect. 0 Within the possible materials, Ti 02 is 
one of the most studied ones. [51. Il5l - [l9j Inside memristive 
devices, T 4 O 7 channels immersed in the Ti 02 -based ma¬ 
trix have been observed, and the formation of such struc¬ 
tures is believed to be largely responsible for the change 
in resistivity. 0 

While Ti 02 is known to be a wide gap semiconductor 
(experimental band gap E g « 3.1 eV[20|) T 4 O 7 presents 
three phases with different electronic transport proper¬ 
ties. The low temperature (T < 140 K) and intermediate 
temperature phases (140 K < T < 150 K) are semicon¬ 
ductors while the high temperature phase (T > 150 K) is 
metallic. mm The differences between those phases is 
mainly due to small atomic displacements which do not 
lead to significant changes in the unit cell. 

Despite the large amount of theoretical work on ei¬ 
ther TiC >2 rutile (l7 L [251 \2l\ and Tup7.[28l-[32j the band 
alignment of these two materials is not known up to now, 
mainly because of the difficulty in building a supercell 
containing the Ti 407 -Ti 02 interface - the two materi¬ 
als present very different unit cells, making the match 
between them very difficult. This information could be 
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very useful to the understanding of the atomic level mech¬ 
anism of the memristor, given that the Ti 407 -Ti 02 in¬ 
terface spans a great part of the device length after the 
formation of the conducting channels. 

As the total energy is given with respect to an arbi¬ 
trary reference energy in periodic calculations, a common 
reference should be determined for all the structures in¬ 
volved in this work.[33| In fact, the knowledge of this 
common reference is necessary to find the offset between 
the valence band of different materials. The knowledge of 
this band offset is important, for instance, in determin¬ 
ing whether there is charge transfer at the interface be¬ 
tween the two systems, the materials doping limits 133 as 
well as other properties. For instance, Scanlon et al. [35[| 
determined the band offset between rutile and anatase 
Ti 02 structures via a cluster approach, providing an un¬ 
derstanding about the separation of charge carriers in 
photocatalysis using the mixed-phase structure. 

Our objective in this paper is to study the interface 
of Ti 02 and one of its oxygen deficient counterparts, 
T 4 O 7 . For this purpose, we built a supercell containing 
the Ti 407 -Ti 02 interface which was then used for our 
calculations. In particular we analyze the band align¬ 
ment, which has been found to be of type II (the conduc¬ 
tion band minimum (CBM) is located in one material 
and valence band maximum (VBM) in the other), and 
its consequences to the overall electronic structure of the 
combined system. 


II. COMPUTATIONAL DETAILS 

All density functional theory (dft[H 0 ) calcu¬ 
lations in this work were performed using the VASP 
package. [HI [39| Wave functions are described within the 
Projector Augmented Waves (PAW) scheme [Id, EU with 
periodic boundary conditions. The 3p3d4s and 2s2p elec¬ 
trons were considered as valence electrons for Ti and O 
atoms respectively and the PBE functional (42j was used 



(a)Ti 4 0 7 (b)Ti0 2 



(c)Ti 4 0 7 -Ti0 2 interface 


Figure 1: (a) Ti 4 0 7 structure generated via the operations given by equations [H and [2j (b) Ti0 2 transformed according to 
equation [T] (c) and Ti 4 0 7 -Ti0 2 interface. The red plane in (c) indicates the interface region and the blue horizontal arrow 
corresponds to the c crystal vector. 


as the exchange-correlation potential. The cutoff energy 
for the plane wave expansion was 700 eV and the k-point 
meshes used throughout this work were generated using 
a 2 x 2 x 1 T-centered grid. Spin polarization was taken 
into account for all calculations. Relaxation of the ions 
was performed while the unit cell was kept fixed, until 
forces were smaller than 2.5 x 10 -3 eV/A, for the iso¬ 
lated materials and 5.0 x 10 -2 eV/A, for the supercell. 
We have performed test calculations with a smaller su¬ 
percell for the interface, but with more strict parameters 
for the convergence of the forces. The results were similar 
for quantities calculated in this work. 

A Hubbard U parameter was introduced for a better 
description of the Ti(d) electrons. Our earlier work [§2] 
shows that a value of 5 eV for U is enough to achieve a 
good description of the electronic structure of Ti n 02 n -i 
(2 < n < 5) Magneli phases. Even though, test calcu¬ 
lations were performed for other values of U (0, 2 and 
4 eV) and the same qualitative results were obtained for 
the type of alignment as well as the relative position of 
the unnocupied levels of the same materials. In the ro- 
tationally invariant scheme as implemented in VASP[43j 
only the value of U — J (J is the spherically averaged 
matrix element of the screened Coulomb interaction be¬ 
tween electrons) is of relevance, thus J was set to zero 
for all calculations. 


III. CONSTRUCTION OF THE INTERFACE 


One of the key points in this work is that building the 
interface requires matching two of the three unit cell vec¬ 
tors of each compound at the interfacial region. For that 
to be accomplished, it is necessary to find the smallest su¬ 
percell where the deformation imposed to each compound 
is such that no significant change in electronic structure 
with respect to the bulk material arises in regions away 
from the interface. For the materials presented in this 
work, this is difficult, since the Ti02 and the Ti 4 0 7 pre¬ 
sented very distinct unit cells. The stable form of Ti02 at 
ambient conditions is the rutile structure (tetragonal unit 
cell, a r = b r = 4.59A, c r = 2.96A[44|) while the Ti 4 0 7 


presents a complicated Magneli structure (triclinic cell, 
see Tabled]). 


Table I: Cell parameters for Ti 4 0 7 , high (HT), intermediate 
(IT) and low temperature (LT) phases [23]. as well as the re¬ 
sult of the transformation given by equation [I] (TR). 



HT 

IT 

LT 

TR 

a (A) 

5.60 

5.60 

5.63 

5.46 

6(A) 

7.12 

7.13 

7.20 

7.14 

c (A) 

20.43 

20.36 

20.26 

20.70 

a(°) 

67.70 

67.76 

67.90 

65.52 

pn 

57.16 

57.44 

57.69 

57.22 

7 n 

108.76 

109.02 

109.68 

108.46 


To overcome this difficulty, we followed steps that 
lead to the derivation of the Magneli structure from the 
rutile one as follows. The first step is to obtain the 
Magneli phase structures from Ti02 rutile via the unit 
cell transformation [45| 
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where a ^, b^ and c^ are the unit cell vectors of the 
Magneli structure of index n. The result of using this 
operation for the case when n = 4 is illustrated in Fig. 
1 1(b) | and the comparison with the Ti 4 0 7 unit cell is pre¬ 
sented in Table [H Notice that this still represents a Ti02 
rutile structure, but the unit cell is now the one from 
Ti 4 0 7 . From this point, the Ti 4 0 7 Magneli phase, as 
well as all Magneli phases (Ti n 02 n -i, n = 2, 3,..., 10) 
can be generated via the operation [ 4614481 ] 

(121) 2- [Oil], (2) 

that consists of periodic displacements; after every n-th 
(121) plane, the atoms are translated by x / 2 [Oil] crystal 
vector of the rutile unit cell. This operation leads to 
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Figure 2: Electrostatic potential profiles along the c vector 
for all structures presented in this work: isolated TiCU and 
Ti 4 Or, and the interface. 



the formation of extended defects, which are composed 
of oxygen vacancies (Vo’s) restricted to planes known 
as shear plan es. The resulting arrangement is presented 
in Fig. |l(a)| Since the operation given by equation [ 2 ] 
affects only the atomic positions, both structures present 
the same triclinic unit cell. 

In essence, the transformation presented above results 
in unit cells with equal lengths. The two materials can 
thus be interfaced i n the plane that is secant to the c axis 
(highlighted in Fig. 1(c)). This Ti 4 0 7 -Ti0 2 interface is, 
by construction, strain free, defect free, presents exactly 
the same cell parameters and consequently no dipoles 
would arise at the heterojunction. In our calculations we 
created a supercell along the c axis by repeating the cell 
for each material four times (c = 165. 6lA) in order to 
ensure convergence of the potential away from the inter¬ 
face. 


IV. BAND OFFSET 


After obtaining the three structures - both isolated ma¬ 
terials and their interface - it is possible to calculate the 
natural band offset between the two systems studied in 
this work. The band offset will be given by |49| 


A^(Ti0 2 |Ti 4 0 7 ) = A£^(Ti 4 0 7 ) - AEy^t (Ti0 2 ) 

+ A£^,(Ti0 2 |Ti 4 0 7 ), (3) 

where AF^(Ti 4 0 7 ) and A£^/(Ti 0 2 ) are the dif¬ 
ferences between the Valence Band Maximum energy 
( Evbm ) and the electrostatic potential averaged over a 
certain volume of the unit cell of the isolated materials, 
while AF^ 5 0 /(TiO 2 |Ti 4 O 7 ) is the difference between the 
averaged electrostatic potential of the different materials 
calculated in bulk-like regions of the relaxed heterostruc¬ 
ture supercell. 

The determination of the quantities needed by equa¬ 
tion [3] can be achieved by, first, obtaining the electronic 
structure of each isolated (bulk) system. From this calcu¬ 
lation, the average electrostatic potential along the crys¬ 
tal axis which is secant to the interface (the c axis, in 


this case) is determined, 

y na rb 

Hz) = ^ J o J o u ( x > y > *) da;d y ( 4 ) 

as well as the position in energy of the VBM, Evbm for 
each system. Using the heterostructure after ionic relax¬ 
ation we also calculated (j){z) in a similar fashion. From 
the plot of the electrostatic potential for all three struc¬ 
tures (shown in figure [ 2 ]) it is possible to define regions, 
labeled Az± and Az 2 , where both the bulk materials and 
the heterostructure present similar characteristics. These 
regions of the heterostructure are considered as bulk-like 
and should have, in principle, the same potential. The 
average electrostatic potential is then obtained by the 
integration of (j)(z) along Az% and Az 2 for Ti 4 0 7 and 
Ti0 2 respectively, and is used to calculate the first and 
second terms of equation [3l From the heterostructure 
calculation, the same quantities are used to obtain the 
final term. 

The third step would comprise the determination of 
the band energy shift due to the hydrostatic pressure. 
In this particular case our heterostructure is completely 
strain free, as it is built by equal unit cells, thus this last 
step presents no contribution. After those steps were 
carried out, the band offset obtained in our case resulted 
in AF^(Ti0 2 |Ti 4 0 7 ) = 2.18 eV. Numerical errors due 
to the choice of a particular V were estimated to be « 
0.3 eV, by comparison with other test calculations per¬ 
formed using other values of U. Even though the value 
of the valence band offset presented a variation with re¬ 
spect to the value of U used, being smaller for smaller C7, 
the overall electronic structure of the interface has not 
changed, and the final result, which is the relative po¬ 
sition of electronic levels of Ti0 2 and Ti 4 0 7 , remained 
practically unchanged. For this reason, we have chosen 
to report only the results for U = 5 eV. 


V. RESULTS 

The projected Density of States (PDOS) was calcu¬ 
lated for the two isolated structures and the results were 
shifted in energy by the band offset. These results are 
shown in fig. [3l From that graph, it is possible to no¬ 
tice that the VBM of the interface is located on Ti 4 0 7 , 
while the CBM is found in the Ti0 2 , resulting in a type 
II heterostructure. 


Given that the CBM is 0.16 eV above the VBM (upper 
panel, Fig. n, it is possible that electrons would leave 
the Ti 4 0 7 structure and arrive at the unoccupied levels 
at the Ti0 2 structure due to the presence of an electric 
field. In this case, the Ti 4 0 7 channels present inside of 
the Ti0 2 matrix that forms the memristor could act as 
donors to the neighboring region of that matrix. This 
specific characteristic does not depend on the choice of 
U. While it is known that the isolated oxygen vacancy 
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Figure 3: PDOS of Ti02 rutile and TLO 7 . Gray shaded areas 
denote the occupied levels (below Ef , which is highlighted 
by the dotted vertical lines) and the two spin components are 
pictured by positive and negative values along the vertical 
axis. The upper panel shows the band alignment between the 
two materials. 
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Figure 4: Projected charge density for the CBM (a) and VBM 
(b) of the Tu 07 -Ti 02 interface after ionic relaxation. The 
blue plane points out the interface region of the structure, 
being composed of Tu 07 to the left and Ti 02 to the right of 
the line. 


Further evidence of the real-space separation of the 
VBM and CBM is given by the projected charge density 
over the frontier orbitals. The projection was performed 
for the last occupied and first unoccupied bands, at the 
T point. As a result, presented in Fi g. [H it is possible to 
notice that while the VBM (Fig |4(b)~| ) is located mainly on 
the Ti atoms of the T 4 O 7 structure close to the interface 
(with a contribution of some of the Ti atoms from Ti 02 ), 
the CBM (Fig. |4(a)| ) is located exclusively on the Ti 02 
segment. 

Some hybridization may be responsible for the contri¬ 
bution coming from the Ti 02 structure to the VBM, as 
both levels present Ti(d) character (see Fig [3]). Another 
explanation for this contribution would be the length of 
the c axis, which is not long enough to describe the sep¬ 
aration of the frontier orbitals. Indeed, calculations per¬ 
formed with a smaller supercell (the supercell c vector 
length was half of the one used in the calculations pre¬ 
sented) resulted in larger charge delocalization. There¬ 
fore, we expect that, in the limit of longer unit cells 
(closer to realistic systems) the charge would be fully 
separated on the two phases. 


VI. CONCLUSIONS 


In conclusion, we have reported the band offset given 
by the average of the electrostatic potential for the 
Ti 407 -Ti 02 heterostructure. The supercell containing 
the junction was built using transformations that re¬ 
sulted in the same unit cell for both materials present 
on this junction. The localization of the VBM and CBM 
was determined to be at the TL 4 O 7 and Ti 02 structures 
respectively, resulting in a type II heterostructure. As the 
band gap is 0.16 eV, our results indicate that the T 4 O 7 
channels inside of the memristor Ti 02 matrix would then 
act as donors for the Ti 02 structure. 


in the Ti 02 rutile acts as a shallow donor, as to our 
knowledge, this same behavior has not been reported for 
the ordered defects present into the Ti 4 C >7 structure. 
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